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1. Introduction 

In QCD it is not the quarks and gluons that are asymptotic states, directly observable in detec- 
tors, but rather colourless composites such as mesons and baryons. This entails the need to describe 
in detail the properties of our final states in terms of their constituent particles, an inherently com- 
plicated task due to the non-perturbative effects of confinement and dynamical chiral symmetry 
breaking. 

The natural framework for this composite description of mesons in the continuum are the 
Schwinger-Dyson and Bethe-Salpeter equations (SDEs and BSEs). Such studies have been exten- 
sively performed in the Rainbow -Ladder approximation which generally reduce to a quenched 
study with non-perturbative effects subsumed into an effective gluonic interaction. While there 
have been several studies attempting to improve upon this simplest of all truncations, in the form of 
including unquenching effects due to quarks, e.g. ||, Q], we here take a different viewpoint. In- 
stead of modelling the contribution from quarks directly, we rather consider the dominant contribu- 
tions arising from the resultant mesonic degrees of freedom. Thus in this talk we focus upon pion- 
cloud contributions [^] to the light meson spectrum in a beyond the rainbow truncation scheme. 

2. Including the Pion-Cloud 

The prescription for including pion degrees of freedom in the SDEs and BSEs in a manner 
consistent with the axial-vector Ward-Takahashi identity (avWTI) was first presented in and 
investigated within the real- value approximation. Further modifications were suggested in [^], 
with the resultant system of equations depicted in Fig. |l|. The first loop diagram of both pictorial 
equations relates to the usual rainbow-ladder, where the infrared suppressed gluon is enhanced 
by the vertex dressing, indicated by YM. To satisfy the avWTI this enhancement is restricted to 
depend on the same momentum as the exchanged gluon; the two dressings are often combined into 
an effective gluon dressing which is consequently modelled. The second diagram that contributes 
to the quark-SDE and meson-BSE represents the back-reaction of the pion onto the quark, i.e. pion- 
cloud effects. This requires knowledge of the quark-pion vertex, which we parameterise in terms 
of the pion Bethe-Salpeter amplitude. Necessarily this couples the two equations (whose form is 
presented in [^, ||]), forming a highly non-trivial system of integral equations. 




Figure 1: The approximated quark-SDE with effective one-gluon and one-pion exchange, together with the 
corresponding BSE. The up-down arrow indicates an averaging procedure of the pion-exchange diagram 
with respect to the dressed/undressed quark-pion vertex. 
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2.1 Yang-Mills Part 

We need to specify the details of our gluon exchange and quark-gluon vertex that model the 
interaction. We choose two different model ansatze that have been employed already in previous 
works. The first model takes fits to the gluon as obtained from numerical solutions of the cor- 
responding Schwinger-Dyson equations. This is taken in combination with results inspired from 
solutions of the quark-gluon vertex for a restricted kinematic section This has the virtue that it 
also provides for generation of a topological charge and hence an anomalous mass contribution to 
the 77, T]' iQ. 

Thus, in choosing a particular kinematic configuration for our vertex dressing we employ a 
rainbow-ladder form of the interaction (thus preserving the avWTI), i.e. 

T™{p,,p2,p,) = y,Z2/%T™{pl), (2.1) 

with quark momenta by p\ and p2 and the gluon momentum 773. Note, however, that ( ^ ) involves 
only the 7v-part of the full tensor structure of the vertex. It has been shown in the analysis of the 



full quark-gluon vertex of Ref. [ ]11| ] that such a model cannot capture all essentials of dynamical 
chiral symmetry breaking. We refer to this as the soft-divergent model interaction. 

The second model that we consider is that of Maris-Tandy In this case, the dressing of 
both the gluon and the vertex is modelled by a phenomenological ansatz that includes the correct 
one-loop UV running and provides dynamical chrial symmetry breaking. We include a study of 
this much used model in the context of pion unquenching for comparison. 

2.2 Pionic Contribtion 

The decomposition of a Bethe-Salpeter vertex function T{p;P)^^'> is well-established in the 



literature, with its form constrained by transformation properties under CPT [|13|]. In particular the 
pion is given by the following form 



Ti{p-P)=Tiys F,{p-P)-ifF2{p;P)-i]/> {p ■ P)F^{p-P) -[f J]F^[p-P) 



(2.2) 



This is of particular relevance for our pionic part of the interaction. We approximate the full pion 
Bethe-Salpeter wave function in the quark-SDE and the kernel of the BSE by the leading amplitude 
in the chiral limit given by 

T{{p-p) = Tiy,^^. (2.3) 

Jn 

Here B^{p^) is the scalar dressing function of the quark propagator in the chiral limit. The effects 
of neglecting the three sub-leading amplitudes have been quantified for a real-value approximation 
in Ref. [^. The great advantage of the approximation ( p3| ) compared to the full back-coupling per- 
formed in Ref. is that we can then fully take into account the quark propagator in the complex 
plane as necessary input into the Bethe-Salpeter equation. Whilst proving to be a simple prescrip- 
tion in itself, it gives rise to the technical challenge of having to evaluate the full normalisation 
condition of the BSE due to the non-trivial momentum dependence of the exchange kernel. We 
discuss this in the next section. 



3 



Pion-cloud effects in the BSE description of mesons 



Richard Williams 



3. Normalisation 



Since we solve homogeneous equations for the Bethe-Salpeter amplitude, their subsequent 
renormalisation comes as an auxiliary condition derived from the inhomogeneous BSE: 



2^^tr 



d'^k 

+ 



3 ( T\{k-Q)S{k + P/l)T{{k,Q)S{k-P/l 
[X%{q, -Q)KfZ"{q, k;P) \xi\ut {k, Q) 



(In) 



where = — is fixed to the on-shell meson mass, the trace is over Dirac matrices and the Bethe- 
Salpeter wave-function x is defined by xit{k;P) = S{k + P/2)ri{k,P)S{k-P/2) . The conjugate 
vertex function F is given by r{p,—P) = CT^ {—p,—P)C^^ , with the charge conjugation matrix 
C = -7274- 



The first term of ( |3.1[ ), independent of the kernel, is easily evaluated since one needs only 
derivatives of the quark propagator. The second term is substantially more complicated due to 
the double integration over the interaction kernel with respect to two four-momenta, k and q. We 
evaluate the integral of ( p| ) and employ finite difference methods to compute the derivative. We 
find that the contribution from the kernel provides important contributions. 



4. Results 

The parameters of for both model interactions were fit to pion observables, with the additional 
constraint of the topological charge for our soft-divergent interaction. In constraining the parameter 
set to meson observables, we find for our model a quark mass of m-^ = 3.4 MeV at jj. = 2 GeV, 
whilst for Maris-Tandy we have = 4.4 MeV. The remaining parameters of the interaction are: 



Soft-Divergent Interaction 


Maris-Tandy Interaction 


di dj dj, ^QCD 
(GeV2) (GeV^) (GeV^) (GeV^) 


CO D nit ^QCD 
(GeV) (GeV2) (GeV) (GeV^) 


1.45 0.1 3.95 0.52 


0.37 1.45 0.5 0.2342 



We calculated a range of meson observables, detailed in Table |^. We observe that the effect of the 
pion back-reaction has only a small impact on the pion mass itself, resulting in a small positive 
or negative shift depending upon the form of the interaction. The impact of including pion-cloud 
effects on the leptonic decay constant is fairly large, with effects of the order of 10%. 

For the remaining heavier mesons, the common trend is that the inclusion of such an unquench- 
ing gives rise to negative mass shifts of 100-200 MeV. Most notable of these are for the rho, where 
we predict that unquenching from the pion-cloud yields a bound-state ~ 100 MeV lighter than in 



the quenched theory, in line with recent lattice simulations [|14|]. 

It is clear, however, that in order to reproduce the rich spectrum of light mesons that we need 
to include spin dependent contributions from the Yang-Mills part of the quark-gluon vertex. It is 
envisaged that this will indeed have a strong impact on the calculated masses of bound-states and 
is the object of future research. 
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Model Employed 




Jji 


Ma 


Mp 


■fp 


Mai 


Ml, I 


M„ 




Maris-Tandy 


w/o pi 
inc. pi 


140 

138+ 


104 

93.2+ 


746 

598 


821 

720 


160 

167 


979 

913 


820 

750 




Our Model 


w/o pi 


125 


102 


638 


795 


159 


941 


879 


493 


949 


inc. pi 


138+ 


93.8+ 


485 


703 


162 


873 


806 


497 


963 


Experiment 


138 


92.4 


400-1200 


776 


156 


1230 


1230 


548 


948 



Table 1: BSE results for a range of mesons in the Maris-Tandy and soft-divergence models employed with 
('incl.') and without ('w/o') the pion back-reaction. Model parameters are tuned such that values marked by 
t are reproduced when the pion-exchange kernel is switched on. Results for rainbow-ladder without pion 
effects are for the same parameter set. 
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